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Abstract 

The evaluation of the nuclear matrix elements (NME) of the two- 
neutrino double beta {2vf3fi) decay and neutrinoless double beta {QvPfi) 
decay using the proton-neutron quasiparticle random-phase approx- 
imation (pnQRPA) is addressed. In particular, the extraction of a 
proper value of the proton-neutron particle-particle interaction param- 
eter, g-pp, of this theory is analyzed in detail. Evidence is shown, that 
it can be misleading to use the experimental half- life of the 2vl3(3 decay 
to extract a value for g^p. Rather, arguments are given in favour of 
using the available data on single beta decay for this purpose. 
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The recent large-scale neutrino-oscillation experiments, Super-Kamioka- 
nde [H, SNO [21, KamLAND 0, CHOOZ [l], have confirmed the existence 
of the neutrino mass. These experiments can only probe the differences of 
the squares of the masses, not the absolute mass scale of the neutrino. On 
the contrary, the neutrinoless double beta decay can probe the abso- 

lute mass scale using the effective neutrino mass, {rriu), extracted from the 
results of the underground double-beta-decay experiments. To extract the 
absolute neutrino masses one needs information about the involved nuclear 
matrix elements El , neutrino mixing , and the associated CP phases 
[H]. As a matter of fact, knowing the underlying nuclear matrix elements 
accurately enough, one can extract from the double-beta experiments infor- 
mation about the CP phases of the neutrino-mixing matrix [Hj. 

One more fundamental piece of information would emerge if the 
decay were detected, namely that the neutrino would be a Majorana particle, 
i.e. an object for whom the particle and antiparticle states coincide. The 
OvPP decay then immediately implies also nonconcervation of the lepton 
number, changing the lepton number by two units. Majorana neutrinos 
are naturally contained in various particle-physics theories going beyond 
the standard model, such as grand-unification theories and supersymmetric 
extensions of the standad model. 

Given the above impressive list of important qualitative and quantita- 
tive neutrino properties, potentially probed by the decay, one can not 
stress enough the importance of a reliable calculation of the involved nuclear 
matrix elements (NME). Lack of accuracy in the values of these matrix ele- 
ments is the source of inaccuracy in the information on the neutrino masses 
and CP phases, extracted from the Oz^/3/3-decay experiments. In particu- 
lar, in view of the planned near-future large-scale underground experiments, 
with detectors in the ton scale, knowledge of the most promising nuclear 
candidates for detection is of paramount importance. 

Contrary to the decay, the two-neutrino double beta {21/(313) decay, 
with two neutrinos and two electrons in the final state, can proceed as a 
perturbative process within the standard model. It can also be used as a 
test bench for the nuclear models, since the decay proceeds via only the 1"*" 
states of the intermediate double-odd nucleus. Success in describing this 
decay mode is a prerequisite for a reliable calculation of the NME's related 
to the Qv(3(3 decay. During the last two decades a host of different nuclear 
models have been used to compute values of the matrix elements involved 
in both types of double-beta-decay transition^ [HI Ej. The mostly used 
nuclear models in the evaluation of the NME's of double beta decay are 
the nuclear shell model and the proton-neutron quasiparticle random-phase 
approximation (pnQRPA), designed for spherical or nearly spherical nuclei. 

After the first shell-model attempts, the problem of the NME's of the 
21/(3(3 and Qv(3(3 decays was viewed in a fresh new way by the introduc- 
tion of the pnQRPA with an adjustable particle-particle part of the proton- 
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neutron two-body interaction. Determination of the value of the correspond- 
ing strength parameter, ^pp, has been a key issue since the mid 80's. As no- 
ticed in the early works jlH I12j , the NME of the 21/(5(5 decay is very sensitive 
to the value of this parameter, leading to the so-called Qpp problem of the 
pnQRPA. On the other hand, the NME of the Qv(5(5 decay is much less de- 
pendent on the value of ^pp, as discussed, e.g., in OEl- Many extensions of 
the pnQRPA have come to light during the last nine years. The first of these 
was the so-called renormalized pnQRPA (pnRQRPA of Ref. ^5)- Other 
extensions of the pnQRPA, used in the /3/3-decay calculations, are cited e.g. 
in 13 El ll'^l llHj . A common feature of all these extensions is the attempt to 
introduce the Pauli exclusion principle into the pnQRPA by improving on 
the quasiboson commutation relations, adopted at the pnQRPA level. In 
these theories different types of correction to the boson commutators of the 
bifermionic operators are introduced, leading to renormalization factors at 
the level of the pnQRPA equations of motion. 

The results of the /3/3-decay calculations are quite scattered ^2] (see also 
[Hlini for a detailed discussion of the matrix elements up to the year 1998), 
and recently it has been suggested that this shortcoming could be overcome 
in the framework of the pnQRPA and its renormalized extensions. In this 
scheme it has been suggested that one could use data on the 2v(5(5 decay 
to extract a more accurate value for the NME corresponding to the neutrino- 
mass mode (i.e., decay mode mediated by the mass of the neutrino) of the 
01/(5(5 decay. The essentials of this method are summarized as follows: the 
value of the interaction strength parameter ^fpp of the pnQRPA (or any of 
its renormalized extensions) can be determined by fitting the value of the 
computed NME to the one extracted from the experimental half-life of the 
corresponding 2i>P(5 transition. This fitted value of gpp is then used in the 
computation of the Qv(5(5 NME. This suggestion has recently been made 
also in In the following, the implications and pitfalls of this scheme 

are analyzed in detail by using the simple and transparent framework of 
the plain pnQRPA. The same qualitative features persist largely also in its 
renormalized extensions. At the same time, arguments are given in favour 
of an other approach, namely fitting g^p by the data on single beta decay (s). 

To have an idea of the suggested procedure JB]) and its alternative, 
advocated in this work, it is instructive to write down an expression for the 
2i//?/3-decay half-life, j for a transition from the initial ground state, OJ^, 
to the final ground state, Op. This expression reads 



a/2 



-1 



(2v) 



where G^'^'^^ is an integral over the phase space of the leptonic variables jSj. 
The nuclear double Gamow-Teller matrix element, M^^q^, corresponding to 



3 



the 2vl3l3 decay, can be written as 



DGT 



(Op iii:,cto>-||i+) 



X 



^ {\Qpp + En- Mi)/me + 1 
(l+||5:^(i)t7l|0+), 



(2) 



where the transition operators are the usual Gamow- Teller operators for (3~ 
transitions, Qpp is the 2v(3(5 Q value, En is the energy of the nth inter- 
mediate state. Ml is the mass energy of the initial nucleus, and rUe is the 
rest-mass of the electron. 

As an alternative to the proposed ^lE] use of the measured 2iyf3(3 decay 
half- life to determine the value of gpp, the use of the measured single-beta- 
decay half-lives is advocated in this work. The available data on Gamow- 
Teller transitions of heavy double-odd nuclei, involved in double and 
double /3''~/EC transitions, have been summarized in Tabled At the mo- 
ment, it is believed that the double /5~ decays are better accessible to ex- 
periments than the double /J'^/EC decays. Nevertheless, it is instructive to 
show the available beta-decay data for nuclei involved in the double 
decays, as well. 

As the first, clean-cut test case one can take the decay of ^^^Cd which 
is a nearly spherical, almost semi-magic nucleus. The corresponding final 
nucleus of the 2i//3/3 decay is ^^^Sn, a genuine spherical semi-magic nucleus. 
Both these nuclei are well describable by the spherical pnQRPA. 

The calculation of the matrix element of Eq. ((2)) proceeds on the following 
lines. The single-particle energies of the sherical mean field are obtained 
from a Woods-Saxon single-particle potential, including the Coulomb and 
spin-orbit parts in the Bohr-Mottelson parametrization j21j . The single- 
particle valence space is taken typically to span two to three oscillator major 
shells around the proton and neutron Fermi surfaces. The adopted two-body 
interaction is a realistic one, based on the one-boson-exchange potential of 
the Bonn type, transformed to nuclear matter by the G-matrix technique. 
The finite-size effects have been taken into account in an approximate way 
by using simple scaling parameters for the short-range monopole part, and 
separate scalings for the = multipole in the particle-hole and particle- 
particle channels. Details of the calculation can be read in j22j . 

The strong short-range correlations between nucleons have been treated 
by using the BCS approximation. The associated pairing strengths are ad- 
justed to reproduce the empirical pairing gaps, extracted from the exper- 
imental separation energies of protons and neutrons, in a way described 
in [231 • The proton-neutron correlations are treated at the pnQRPA level 
by fixing the scale of the particle-hole J'^ = 1^ two-body matrix elements 
to reproduce the empirical location of the Gamow-Teller giant resonance, 
whereas the particle-particle part of the same interaction is scaled by the in- 
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Table 1: Experimental EC- and (3^ — decay log ft values for heavy double- 
odd nuclei involved as intermediate nuclei in double /3~ and double 
decays. For completeness, also the Q values of the /?/? decays are given in 
the second column. 



(3(3 mode 


Q [MeV] 


Init. nucl. 


Final nucl. 


Mode 


log/t 


Ref. 


rr 


3.03 


looxc 


iooMo 


EC 


4.45 








iootc 


iooRu 




4.6 


d 


rr 


1.30 






EC 


4.3 


d 








104p^ 




4.5 


in 


/?+/?+ 


0.73 


106 Ag 


106p^ 


EC 


4.9 


IE] 






106 Ag 


lo^Cd 




> 4.2 




p-(3- 


2.01 


iiOAg 


llOpd 


EC 


4.1 


m 






iiOAg 


iiocd 




4.7 


nu 


(3-(3- 


0.53 


114j^ 


ii^cd 


EC 


4.9 


uni 






114j^ 


ii^Sn 




4.5 


nni 


rr 


2.80 




ii^Cd 


EC 


4.39 


uni 












4.7 


m 




0.87 


128j 


i28Te 


EC 


5.0 


m 






128j 


i28Xe 




6.1 


m 




0.54 


130Cs 


i30Xe 


EC 


5.1 


m 






130Cs 


i30Ba 


/?- 


5.1 


m 




0.37 


i36La 


i36Ba 


EC 


4.6 


m 






i36La 


136 Ce 


/3- 


? 


m 



teraction strength constant (/pp, left as a free parameter in the calculations. 
This method was used for realistic interactions in the context of the 2v(3(3 
decay in j^, and in description of single beta decays in pS] . 

In Fig. HI panel (a), the NME M(2'^)(tot), corresponding to the 2v(3(5 
decay of ^^^Cd, is drawn as a function of ^pp. In the same figure a rough 
value of the extracted experimental NME, M(^'^)(exp.), has been shown as 
a horizontal line, since its value is independent of ^pp. Here the uncer- 
tainties in the value of this extracted NME, arising from the experimental 
error in the measured half-life, and the uncertainty in the proper value of 
the axial-vector coupling constant, (ja, for medium-heavy and heavy nuclei, 
have been omitted. The intersection point of these two curves gives now the 
fitted value, (7pp(/3/3) ~ 1.03, of g^p. As can be seen from the curve denoted 
by M(^'^^(1^) in Fig. ^ the NME including only the contribution arising 
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Figure 1: Panel (a): The NME's corresponding to the 2v(5(3 decay of ^^^Cd 
shown as functions of ffpp. The complete NME, M(2'^)(tot), the NME with 
only the lowest intermediate contribution included, M^'^'^^V^), and the ex- 
perimental NME, M(^'^''(exp.), have been shown. Panel (b): The left-branch, 
EC NME, Mec, and the right-branch NME, /?" NME, Mp-, shown as func- 
tions of g'pp. Panels (c) and (d): The same as (a) and (b) for the 2vl3l3 decay 
of i28Te. 

from the virtual transition through the first l"*" state, 1^, of the interme- 
diate nucleus -"^-^^In, almost coincides with the complete NME, M(^'')(tot), 
especially for g'pp values around unity. This is a characteristic of the so-called 
single-state dominance (SSD), studied extensively e.g. in [21]. 

In the case of such a SSD, the NME ^ of the 2i>(5(5 decay can be 
approximately written as 

M(-) _ MjeMs_^ . (3) 

(iQjj + El - Afi)/m.o + 1 

The two branches of the 2i^(3(3 transition, Mec and Mg- are drawn as func- 
tions of Qpp in panel (b) of Fig. ^ It is remarkable that the magnitudes 
of the left-branch NME, corresponding to the electron-capture (EC) decay 
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of the 1^ state in ^^^In to the ground state of ^^^Cd, and the right-branch 
NME, corresponding to the j3~ decay of the same state to the ground state 
of ^^^Sn, can in some cases be determined from experimental data on the 
corresponding decay half-hves. It is also clear that in this kind of a simple 
case the study of the relation between the single and double beta decays 
is most transparent, in particular, related to the determination of the g^^ 
parameter. 

Using the extracted value of gpp{(3(3), one immediately obtains, due to 
the SSD, the values of the left- and right-branch NME's, as shown in panel 
(b) of Fig. HI From Fig. [Hone obtains Mec ~ 1-4 and M^- ~ 0.24. These 
values of the NME's can, in turn, be used to compute the half-lives of the EC 
and P~ decays from the 1^ state in ^^^In. Comparison of these computed 
values with the corresponding experimental ones, extracted from Table ^ 
yields 

4^,^)(exp.) t[%\exp.) 

^PP^ ^ 2.6 ; ~ 0.16 , (4) 

4^C)(th.) t[%\th.) 

indicating that for gpp{(3(3) ~ 1.03 one obtains too fast an EC transition 
and much too slow a (3~ transition. Fitting the f3~ decay half- life, instead 
of the 21/(3(3 decay half-life, would yield a value gpp{(3~) ~ 0.85, which also 
would result in a more reasonable matrix element for the EC branch, namely 
Mec — 1-2. The corresponding experimental magnitude is MEc(exp) ~ 
0.8, the exact value depending on the adopted value for g^. As can be 
seen, the proper determination of the value of the i^pp parameter, by using 
the (3~ decay half-life, can lead to a notably different value from the one 
extracted by using the 2^(3(3 decay half-life, even in the simple case of the 
SSD. Summarizing the above: use of the value gpp{(3(3) ~ 1.03 reproduces 
the 2v(3(3 half-life via two compensating errors: too large an EC NME is 
compensated by too small a (3~ NME. 

As the second test case one can take the 2u(3(3 decay of ^^^Te to the 
ground state of ^^^Xe. This case can be analyzed using the very methods 
deviced for ^^^Cd in Fig. ^ panels (a) and (b). A corresponding scheme is 
shown for the ^^'^Te decay in Fig. Q panels (c) and (d). As can be seen from 
panel (c), the curves for the total matrix element and the M(^^)(1^) matrix 
element are very much separated everywhere but at the values of gpp close 
to the point which reproduces the value of the experimental matrix element. 
Hence, in this case one can not speak about SSD, and the situation is more 
complicated than in the ^^^Cd case. In this case the intersection point of the 
curves, corresponding to the total and experimental matrix elements, gives 
the fitted value, gpp{(3(3) ~ 0.82, of gpp. 

The two branches of the matrix element M^"^^") (V^) , Mec and Mg-, 
are drawn as functions of ^pp in panel (d) of Fig. ^ Using the extracted 
value of gpp {(3(3) , one obtains for the left- and right-branch matrix elements 
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Mec — 1-19 and Afg- ~ 0.05. These values of the NME's can, in turn, 
be used to compute the half-hves of the EC and (3~ decays from the 1^ 
state of ^^^I. Comparison of these computed values with the corresponding 
experimental ones, extracted from Tabled yields 

t'f,'^'(exp.) _ t'f-\e^p.) _ 

indicating that for gpp(/3/5) ~ 0.82 one obtains much too fast an EC tran- 
sition and much too slow a (3~ transition. Fitting the f3~ decay half-life, 
instead of the 2i//3/3 decay half-life, would yield a value gpp{(3~) ^ 0.755, 
which would only slightly change the value of the matrix element for the 
EC branch, namely to Mec ~ 1.15. The corresponding experimental mag- 
nitude is MEc(exp) ~ 0.38, for = 1-0. As can be seen, in this case the 
proper determination of the value of the g^p parameter, by using the (3~ 
decay half-life, does not lead to a notably different value of Mec from the 
one extracted by using the 2z//3/3 decay half-life. The reason for this discrep- 
ancy is not clear, but deformation effects could play some role. Even so, the 
above tells us that the use of the value QppiPP) — 0.82 reproduces the 2i//3/3 
half- life via two compensating errors: too large an EC NME is compensated 
by too small a /3" NME. 

As the third case, the 2i//3/3 decay of ^^Ge to the ground state of ^^Se 
will be discussed. This case can be analyzed along the lines of the previous 
two cases. A corresponding scheme is shown for the '^^Ge decay in Fig. |2l 
panels (a) and (b). As can be seen from Fig. [2 there exists no SSD, and 
the situation is in this respect similar to the ^■^^Te case. In this case the 
intersection point of the curves, corresponding to the total and experimental 
matrix elements, gives the fitted value, gpp{f3f3) ~ 0.94, of gpp. 

The two branches of the matrix element M^^'^''(l^), Mec and Mg-, are 
drawn as functions of gpp in panel (b) of Fig. [3 Using the extracted value 
of gpp{(3f3), one obtains the values Mec — 1-52 and Mg- ~ 0.09 of the left- 
and right-branch NME's, which give for the corresponding log ft values 

log/t(EC)~3.9 ; log/t(/3-) ^6.4 , (6) 

for fl'A = 1-0. The lowest state in the intermediate nucleus, ''^As, is a 2~ 
state, and hence the Gamow-Teller decays of the lowest l"*" state are hard 
to observe due to the fast gamma decays to this 2~ state. 

As the next example of the gpp{/3/3) problem, the 2vf5[5 decay of *^Se to 
the ground state of ^^Kr is discussed in Fig. [21 panels (c) and (d). As for 
the ^^Ge case, also here the SSD is not applicable. From Fig. [Jlone can read 
for the intersection point of the total and experimental matrix elements 
the value 5pp (/?/?) — 1.07, giving for the EC and /3~ NME's the values 
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Figure 2: The same as Fig. ^ for NME's corresponding to the 2v(3(3 decays 
of ^^Ge and ^^Se. 

Mec — 1-32 and Mp- ~ 0.11. These, in turn, give for the corresponding 
log ft values 

log/t(EC)~4.0 ; log/t(/3-) -6.2 , (7) 

for = 1-0. The lowest two states in the intermediate nucleus, ^^Br, are a 
5~ state and a 2~ state, and hence the Gamow-Teller decays of the lowest 
l"*" state have not been observed. 

Although no measured EC or (3~ NME can be extracted for the ^^Ge and 
^^Se cases, one can compare the computed log ft values of Eqs. © and Q 
to the log ft values of similar cases in the same mass region. In the relevant 
mass region there are three double-odd nuclei with a l"*" ground state and 
decay patterns analogous to the ones of ''^As and '^^Br, namely the ones 
listed in Table 121 From this table one immediately notices that the log ft 
values of the relevant EC decays range between log /t(EC) = 4.7 — 4.8 and 
values of the relevant f3~ decays range between log ft{f)~) = 5.1 — 5.5. This 
would suggest that the extracted log ft values of Eqs. © and for the 
EC decays are too small and the corresponding extracted log/i values for 
the f3~ decays far too large. Much better agreement between the theoretical 
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Table 2: EC- and /5 — decay log ft values for selected decays of double-odd 
nuclei in the pf shell. The data is taken from |19j . 



Init. nucl. 


Final nucl. 


Mode 


log/t 


70Ga 


™Zn 


EC 


4.7 


^OQa 


™Ge 


/3- 


5.1 


78b, 


78Se 


EC 


4.8 


78Br 


78Kr 




? 


80Br 


80Se 


EC 


4.7 


80Br 


soKr 


P- 


5.5 



and experimental EC and f3^ log ft values, around log /t(EC) ~ 4.6 and 
log/i(/3~) ~ 5.3 could be obtained for smaller values than the one, 
suggested by the 2i//3/3-decay half-life. For the two discussed decays a value 
5pp(/3~) — 0.8 would do quite well. 

Based on the previous analysis one can say that the conclusions arising 
from the analysis of the 2uf3f3 decays of ''^Ge and ^^Se coincide with the ones 
arising from the 21/(3(3 decays of ^^^Cd and ^^^Te: cancellation of errors in 
the two matrix elements, Mec and Mp-, conspire to produce a 2v(3(3 NME 
which exactly reproduces the corresponding experimental matrix element. 
This demonstrates that it can be dangerous to determine the value of g^^ 
by fitting the 2v(3(3 decay half-life. 

In fact, determination of the value of g^p by the data on single beta 
decay leaves the 2z//3/3-decay half-life as a prediction of the theory. Com- 
parison of this prediction to the experimental half-life would tell about the 
predictive power of the adopted theoretical framework, in terms of the size of 
the adopted single-particle space, the adopted single-particle energies, etc. 
A rougly correct prediction for the 2z//3/3-decay half-life would shed more 
confidence on the theoretical predictions concerning the other multipoles, 
involved in the Qv(3(3 decay. 

As the final example of the gp^{(3(3) problem, the 2v(3(3 decay of -'^'^^Mo 
to the ground state of ^'^'^Ru is discussed in Fig. |3J In this case the SSD is 
roughly applicable. From Fig. ISlone can see that in this particular nuclear- 
structure calculation, the one of Ref. |^ where one can read more details of 
the used single-particle basis, etc., the computed total NME never reaches 
the experimental NME, extracted by using gx = 1.0. Hence, in this case one 
is forced to use the experimental /3~-decay log ft value, quoted in Table Q 
to determine the value of g'pp, resulting in gpp{(3~) ~ 1.02. This gives for 
the EC the value Mec — 1-97, and for the corresponding log ft value 

log/t(EC) -3.7 , (8) 
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0.0 H > 1 > 1 > 1 > 1 > 1 

0.5 0.6 0.7 0.8 0.9 1.0 1.1 




0.0 H • 1 • 1 • 1 • 1 • 1 

0.5 0.6 0.7 0.8 0.9 1.0 1.1 
gpp 



Figure 3: The same as Fig. ^ for NME's corresponding to the 2v[3l3 decay 
of i°°Mo. 

using gp^ = 1.0. This is too low a value for this decay, as seen from the data 
of Tabled indicating that some nuclear-structure effects, e.g. deformation, 
beyond the reach of the spherical pnQRPA, might be present. 

Summarizing the above presented results, the problem of determination 
of the proton-neutron interaction strength, (7pp, in a pnQRPA type of cal- 
culation, be it the plain pnQRPA or one of its renormalized extensions, has 
been addressed. The apparent solution of the "g'pp problem" by fitting i^pp 
to available data on 2i^/3/3-decay half-lives has been critically analyzed. Fit- 
ting this parameter to the existing data on single (5~ transitions is found to 
be a more meaningful solution to the problem. Arguments favouring this 
method have been summarized in TableElwhere the positive points (-I-) and 
negative points (— ) of the two fitting methods have been listed. Below few 
comments concerning the listed points of the table are made. 

Concerning point one, one can even perform a systematical study of the 
beta-decay properties of a given nuclear region in the fit to beta decays. 
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Table 3: Pros (+) and cons (— ) of the two discussed recipes to fit the 
parameter gpp. For more explanation on the various points see the text. 



Point 


Fit to (3" and/or EC decay (s) 


Fit to 2^(3" (3' decay 


1 


One, two or more observables 


Only one observable can be 




can be used for the fit (+) 


used for the fit (— ) 


2 


Direct access to grass-root-level 


Two or more compensating 




deficiencies of a nuclear model (+) 


errors may conspire to produce 
a good 2vj3~ I3~ decay rate (— ) 


3 


The beta-decay properties better 


The 2vl3~ I3~ decay properties 




reproduced (+) 


better reproduced (-I-) 


4 


Error limits from comparison 


Advisable to check against 




of the experimental and computed 


data on j3~ decays 




2i'/5~/3~ decay rate (+) 




5 


Largely eliminates the model-space 


Largely eliminates the model-space 




dependence of the computed 


dependence of the computed 




Ou(3~P~ decay rates (+) 


Oi'P~P~ decay rates (+) 


6 


Can be extended to study of 


No access to a possible variation 




forbidden contributions. 


of gpp from multipole 




e.g. 2~, in OuP~P~ decay (-I-) 


to multipole (— ) 


7 


Can access decays where no 


Can access /3/3 decays where no 




2u(3(3 data exists, see 


direct /3-decay data exists 




Table 121 (+) 


(™Ge and ^^Se) (+) 


Balance 


7x(+) 


3 X {+) and 3 x (-) 



This approach would correspond to a shell-model _24] type of application 
of the beta-decay data. Referring to point two, a fit to (3~ data can reveal 
deficiencies in the predictive power of the used nuclear model in the case of 
the EC rates of the other branch. This seems to be the case, e.g., in the 
present calculation. 

In regard to point four, the first method can be used to draw some 
conclusions about the error limits in the (3(3 calculations, whereas in the 
second method one necessarily should check the consistency of the calcu- 
lations against the available beta-decay observables. This is a necessary 
procedure, not warranting either a plus or a minus mark. In point five the 
similar behaviour of the two discussed fitting methods comes, on one hand. 
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from the fact that in both methods one fixes first the pairing parameters by 
semiempirical pairing gaps. This is an essential step and produces, for each 
single-particle space, the consistent quasiparticle mean field. On the other 
hand, as the next step, both methods use experimental data to fit the i^pp 
parameter. This two-step fitting procedure is enough to eliminate almost 
completely the dependence of the computed Oi^/5/5-decay rates on the size of 
the model space. 

The point number six is a very important one considering the actual 
computation of the Oz^/3/3-decay rates. In the first procedure a separate ^pp 
analysis of higher multipoles can be performed, e.g., in the pf shell where 
data on beta decays of 2~ states are available. In the second method the 
same value of the gpp parameter has to be assumed for all multipoles. Finally, 
it is to be noted that in the previous analysis the axial-vector coupling 
constant, gA, has been assumed to be roughly the same for both the (3 and 
f3f3 decays. Since no exhaustive studies of this matter have been performed, 
we take this assumption at face value in this work. 

Concluding, the last line of Table 01 sums up the positive and negative 
points of each method. This final balance clearly supports the argument 
that the beta-decay fitting should be favoured, rather than the 2z^/?/?-decay 
fitting. 
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